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摘 要: 主 元 分 析 (PCA) 与 二 维 主 元 分 析 (2DPCA) 是 两 种 典型 的 图 像 特征 提取 方法 ， 它 们 所 提取 的 图 像 
特征 的 优 劣 可 由 重建 误差 来 衡量 。 通 过 对 PCA 和 2DPCA 的 重建 误差 分 析 发 现 ， 二 者 的 重建 误 
差 在 理论 上 相同 ， 但 在 实际 应 用 中 取决 于 它们 的 样本 协 方差 阵 的 估计 准确 度 。 本 文 以 均 方 误差 为 
度量 给 出 了 PCA 5 2DPCA 样本 协 方差 阵 的 估计 准确 度 表 达 式 ， 并 由 此 得 到 2DPCA 图 像 特征 
优 于 PCA 的 判定 条 件 是 2DPCA 协 方差 阵 的 特征 值 平方 和 大 于 PCA。 本 文 还 指出 行 2DPCA 与 
列 2DPCA 所 提取 的 图 像 特 征 熟 优 熟 劣 也 取决 于 它们 各 自 协 方差 阵 的 特征 值 平方 和 的 大 小 。 在 人 
脸 图 像 库 与 人 脸 表 情 图 像 库 上 的 实验 验证 了 上 述 判 定 条 件 的 正确 性 。 
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主 元 分 析 (Principal Component Analysis, PCA) 与 二 维 主 元 分 析 (2-Dimension PCA, 2DPC 
A) 是 两 种 典型 的 图 像 特征 提取 方法 。PCA 最 早 于 1901 Eid t DI, JT 20 世纪 90 年 代 前 后 
引入 到 人 脸 图 像 的 表示 与 识别 中 3; 2DPCA 最 早 于 2004 年 提出 并 用 于 人 脸 图 像 表 示 与 识 
别 向 。 以 下 将 由 PCA 从 图 像 中 提取 的 特征 称 为 PCA 图 像 特征 ， 由 2DPCA 从 图 像 中 提取 的 特 
征 称 为 2DPCA 图 像 特 征 。 

PCA 与 2DPCA 的 原理 本 质 上 一 致 ， 区 别 在 于 PCA 需要 先 将 图 像 拉 直 成 向 量 再 处 理 ， 
而 2DPCA 不 需要 这 个 过 程 直 接 对 二 维 图 像 进行 处 理 。 在 进行 图 像 特征 提取 时 ，PCA 5 2DPCA 
都 需要 处 理 一 个 协 方差 阵 来 获取 投影 矩阵 。 上 述 的 区 别 导 致 它们 要 处 理 的 协 方差 阵 有 很 大 不 
同 。 对 于 一 个 m x n 的 图 像 ，PCA 要 处 理 的 是 一 个 rm 阶 协 方 差 阵 ， 称 之 为 PCA 协 方差 阵 ; 
而 2DPCA 要 处 理 的 是 一 个 m (或 m) 阶 图 像 协 方差 阵 咎 ， 称 之 为 2DPCA 协 方差 阵 。 在 实际 问 
题 中 ，PCA 与 2DPCA 的 协 方差 阵 都 由 它们 相应 的 样本 协 方差 阵 来 估计 。 

2DPCA 协 方差 阵 实 际 上 是 由 PCA 协 方差 阵 的 部 分 元 素 构成 的 ， 因 而 它 含有 的 信息 量 不 
如 PCA 协 方差 阵 多 。 但 是 实验 表明 外 ，2DPCA 图 像 特 征 在 图 像 表 示 与 图 像 识 别 两 方面 都 优 
于 PCA 图 像 特征 。 本 文 试 图 对 此 给 出 合理 解释 。 首 先 对 PCA 5 2DPCA 图 像 特征 的 重建 误差 
进行 了 分 析 ， 结 果 表 明 它 们 的 重建 误差 在 理论 上 是 相同 的 ， 但 在 实际 中 取决 于 各 自 样本 协 方差 
阵 的 估计 准确 度 。 本 文 以 均 方 误差 为 度量 对 PCA 与 2DPCA 样本 协 方差 阵 的 估计 准确 度 进行 了 
分 析 ， 得 到 了 2DPCA 样本 协 方差 阵 的 估计 准确 度 优 于 PCA 的 一 个 判定 条 件 。 由 类 似 方法 还 可 
得 到 C2DPCA (Column 2DPCA) 5& R2DPCA (Row 2DPCA) RRAZ KHER. 


收 稿 日 期 : 2008-03-10. 作者 简介 : ER (1972 年 11 月 生 )， 男 ， 博 士 . 研究 方向 : 信号 、 信 息 处 理 与 图 像 处 理 . 
* 基 金 项 目 : 国家 自然 科学 基金 (60872084); 高 等 学 校 博 士 学 科 点 专项 科研 基金 (SRFDP-20060003102). 
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2 PCA 52DPCA 图 像 特征 及 其 重建 误差 


PCA 5 2DPCA 图 像 特 征 主要 用 于 图 像 表 示 与 图 像 识 别 ， 图 像 表 示 的 效果 常用 重建 误差 来 
描述 ， 而 图 像 识 别 率 也 与 重建 误差 有 密切 关系 。 因 此 ， 先 对 PCA 与 2DPCA 图 像 特征 的 重建 误 
差 进 行 分 析 。 

设 入 个 m x n BRIERE Ai, 42, :… , AN 是 来 自 于 总 体 4 的 独立 同 分 布 随机 图 像样 本 ， 它 
们 的 列 拉 直 向 量 为 al, a2, , aw Mao 

2.1 PCA 协 方差 阵 与 图 像 特 征 

随机 图 像 矩 阵 ( 即 总 体 ) 4 或 随机 图 像 向 量 a 的 协 方差 阵 定义 为 回 ， 


E = (0i)mnxmn = E[(a — E(a))(a — 已 (ao)) 门 ， (1) 


'"CAÉPCAWZEER, EE TERA ARRA 83 ZRIBABASE. WX ETUR E E, 
特征 值 为 Xi > A. o > Amn > 0， 相 应 的 单位 正 交 特征 向 量 为 pl, po, … , pmn。 取 
前 上 (1 < k < mn) 个 特征 向 量 组 成 矩阵 P = (p pr- ,pp)， 并 令 b = PT(a — E(a)), Hi 
F E(bbT) = diag( 和 A1, Ào, --- , 和)， 包 含 了 a 的 协 方差 阵 的 前 个 最 大 特征 值 ， 故 称 5 为 随机 向 
量 a 的 k 阶 主 元 分 量 昌 。b 就 是 PCA 图 像 特征 。 

在 实际 问题 中 ， 习 由 其 样本 协 方差 阵 


"E " zd 
$= sl Sousa | 


进行 估计 ， 其 中 元 = 三 Xa XI ÈRES inBU kA RCACRS AES PD SL US 8r IE ADAE 
EAREN Ê = (Pi pa, ,次 )。 

2.2 2DPCA 协 方差 阵 与 图 像 特征 

2DPCA 有 两 种 形式 ， 一 种 称 为 列 2DPCA (Column 2DPCA, C2DPCA)， 另 一 种 称 为 行 2D- 
PCA (Row 2DPCA, R2DPCA). 

C2DPCA 97; El sg XA, 


T. = E[(A - E(A)(4 - E(A)] = J” cov(ai), (2) 


其 中 al a2, … ,an 是 4 的 列 向 量 。T。 是 4 的 各 个 列 向 量 的 协 方差 阵 的 和 ; 
R2DPCA 协 方 差 阵 定义 为 向 ， 


T, = E[((A— E(A)) (A - E(A))] = 2, cov(Bi), (3) 


其 中 B1, 05, ---, Bm 是 4 的 行 向 量 。T; 是 4 的 各 个 行 向 量 的 协 方差 阵 的 和 。 

C2DPCA 协 方差 阵 与 R2DPCA 协 方 差 阵 统称 为 2DPCA 协 方差 阵 。 不 难看 出 ，C2DPCA Hh 
Ji: EE Tz 只 保留 了 4 的 同 列 像素 之 间 的 相关 信息 ， 丢 弃 了 不 同 列 像 素 之 间 的 相关 信息 ; R2DP 
CA 协 方差 阵 T; 只 保留 了 4 的 同行 像素 之 闻 的 相关 信息 ， 丢 弃 了 不 同行 像素 之 间 的 相关 信 
息 ， 因 此 工 。 与 T, 都 是 由 学 的 部 分 元 素 构成 的 。 以 下 仅 对 C2DPCA 进行 讨论 ， 所 得 结果 缘 适 
用 于 R2DPCA。 
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对 C2DPCA 协 方差 阵 工 。 作 特征 值 分 解 ， 特 征 值 为 4 > 72 > … > mm > 0， 相 应 的 单位 正 
交 特 征 向 量 为 g, q2, … , qm。 取 前 1 (1 < 1 < m) 个 特征 向 量 组 成 矩阵 @ = (Q, qz, … , qD 
Jy REPE A 的 1 阶 主 元 分 量 为 B= QT7(4 一 EE(A))， 它 就 是 C2DPCA 图 像 特征 。 

在 实际 问题 中 ，T。 由 其 样本 协 方差 阵 


N 
A 1 2 AT 
ng NE 


Hiit, A= h Y Ais XI P. MERED E, WIRI LA BCKBHIERUPDSERE R RAE LE SURE E 


向 量 矩 阵 为 @ = (di do. d)- 
2.85 图像 特 征 的 重建 误差 


先 给 出 三 个 引 理 : 
SI 设 z 为 n 维 随机 向 量 ，E(z) = pu, cov(r) — X. H Anit, Ru 
E(zT Hx) = tr(HY)+ p Hy, (4) 
其 中 tr(.) 表示 矩阵 的 迹 。 
引 理 2 四 ik Bxk，Ckxt，Dixn 为 三 个 矩阵 ， 则 
Vec(BC D) = (DT & B). Vec(C), (5) 


HH SXI; Kronecker 乘积 ，Vec(.) 表示 和 矩阵 的 列 拉 直 运 算 。 
3]283/7 i Dn WIKRE, U 为 其 前 大 个 最 大 特征 值 所 对 应 的 单位 正 交 特 征 向 量 矩 阵 ， 
则 对 于 任意 n x 天 列 正 交 矩阵 站， 有 


tr(U7 DU) = max tr(V7 DV). (6) 


设 由 PCA 图 像 特征 6 对 随机 图 像 向 量 a ETT RE RE 3A SUV a = E(a)-)-Pb, B C2DPCA Kl 
像 特征 BB 对 随机 图 像 矩 阵 4 进行 重建 的 表达 式 为 4 = E(A) + QB, WE 


定理 1 当 l 
2^ = » 
M, E(IA — A|?) = Etlla — al?)- 
证 明 由 引 理 1，PCA 图 像 特征 2 的 重建 误差 为 
E(la — al?) = E((a — E(a))" (1 — PP")(a — E(a))] = tr[(1 ~ PP7)cov(a)] 


= tr(Z) - tr(PTEP) = tr(5) e à Y Ai / Y X). (T) 


{=i 
由 引 理 1 和 引 理 2，C2DPCA 图 像 特 征 B 的 重建 误差 为 
E(|A — Al?) = E([Vec(I — QQT)(A — E(A))]T[Vec(I5, — QQ7)(A — E(A))]) 
= E{(a — E(a))? (I, & (Im — QQ"))(a z E(a))} 


= tr[In & (Im — QQT)cov(a)] = tr(£) — tr((14 8 Q)TX(, & Q)) 


tr(X) — tr(QTT.Q) = tr(Y) (1 一 $: op 30 (8) 
i=1 一 


zl 
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HFT.) =t) MEA 实 天 ， 所 以 结论 成 立 。 
i—1l i=l 
称 


k mn l m 
3X YA 或 2 n 3m 


为 投影 特征 值 比例 。 定 理 1 说 明 ， 在 理论 上 ， 只 要 投影 特征 值 比例 相同 ，PCA 图 像 特征 
5j C2DPCA 图 像 特征 就 具有 相同 的 重建 误差 ， 因 而 它们 的 图 像 表 示 效 果 与 图 像 识别 率 也 大 体 
HE. 

但 在 实际 问题 中 ， 重 建 表达 式 (7) 5 (8) 中 的 已 @ 分 别 由 已 , @ 蔡 代 ， 由 类 似 于 式 (7) 与 
式 (8) 的 推导 过 程 ， 可 得 实际 问题 中 的 重建 误差 为 


E(|la — à|?|P) = tr(X) — tr(ÊTEÊ), (9) 
E(|A — Aj?]Q) = tr() — tr(QTT.Q). (10) 


由 式 (9) 55x (10) 及 引 理 3 可 见 ， 在 实际 问题 中 PCA 5; C2DPCA 图 像 特征 的 重建 误差 取决 
T PCA 样本 协 方差 阵 立 估计 了 的 准确 度 和 C2DPCA 样本 协 方差 阵 工 。 估计 工 。 的 准确 度 。 


3 ”样本 协 方 差 阵 的 准确 度 分 析 


采用 样本 估计 的 均 方 误差 作为 描述 PCA 样 本 协 方差 阵 萝 估计 于 和 C2DPCA 样本 协 方差 
Wc P. fiit T. 的 准确 度 度量 。 

31 样本 协 方差 阵 的 均 方 误差 

为 叙述 的 方便 ， 记 随机 图 像 矩 阵 4 的 第 (站 元 为 必 ， 已 (al) = ul. (4,0) 765 (j, k) TR 
方差 为 of = E((od — uj)(o& — ut). (i1, Ih) 元 、(iz, l2) 元 、(ia, 13) 元 和 (i4, la) 元 的 混合 四 阶 
中 心 矩 为 

Vins. = E[(oi — ui) (o — ui) (ois — ios = mi): 

不 难看 出 ， 协 方差 ott EMIA X H (i+ m(L— 1), j - m(k — 1) 76. 

定理 2 POA FESHI AE X B5 5 RED 


(—-N5 --AN* — 11N? +15N — 6) C cs c 
+ S > > Te a) 


HF |- |r 表示 和 矩阵 的 Frobenius 范 数 。 
证 明 记 样 本 A, FOR (ij) TA al, bi, = ol, aj, Kb 


N 
EES 812,4, 
1 Nx t 
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JW 7; 25 olk 的 估计 为 
由 于 
有 


EDE = E(A ad, E 


_ (N-11 
= C Eoi 


+ 一 1)? 
+> Eloi is 


I 
alk 5 l l pk 
Tij ac N bisbis, 
Sl 


EQ.) = Tal, 
1 N-—1 k k 1 a k k 2 
Hi Le (ol, 一 Hi 208 (= lo 2s uj) m N (ars; s Hj ) 
MT 9155 
1 X 2 
= ui) (o5, EN uy ] TP ATE NzE( 2 (ol, 一 Hi ib ( T (ay。 一 s )) 
-Pr 4424 
N 0? N 
iA ( Y (ob, ib) + a obs, - P (Yo (ol, aD) 
s2=1 s1=1 
S2 天 s S1 天 s 


X KN 7 p(at, — dab - aD 32. - 8) ( 32 6b. - ub) 


(N - 1) Il y 
E poen 


31-71 s2=1 
S1 天 s s2 天 s 
2(N —1P n kk , 4(N— 1? 
Ne NND cio Ojj + ( N4 (cij) 


X E( Y; Y; (ol. — aola. - D) ( 32. V2 (ob — Doa — ub) 


s1=1 t1=1 


s1 天 s t1 天 s 
N? — AN? - 6N — 


N3 


NN 
E( 3 Y utu) = 2 
3 一 1 


i 


= 


c 


32=1 t2 一 1 
s2 天 s t2 天 s 
3 N-1(2N-3 2(N 2N 
vir CODEN D olog + 2) 


DIO E) 


3 一 1 t—1 s1=1 
tzs 81s 
N 
N-—1,x k 1 k k 
-p Cis- ri) N 2.05. u5)) 
527s 


(A— (oi. 一 nx (ain 一 ui) ) 


NE 
txt 


«(Nghe - b - 3 $c. b) 


1221 
tozt 


 .(N- 1)2N - 3) ite 二 (N — )(-N3 -5N? — 8N +6) NN 


N? 2323 N4 iiO jj 


" (N — LD(N5 — 2N* — 2N? + 8N? — 9N + 6) 


Na (o^. (13) 
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因此 


l 
M: 
WE 
M 
M 
E 
E 
S d 
a 
S&F 
t 


E|* — xz 
- : HE TEJE (Et? - aatia + oth)’) 
E ( 3305 十 > » bf, MU) 
1 s—1 


将 式 (12) fl (13) 代入 上 式 再 整理 即 证 。 
定理 3 (i) C2DPCA 样本 协 方差 阵 工 。 的 均 方 误差 为 


2m m n n 
gir. rJ = CS Y Y Yr 
一 N4 十 3N3 NE ANM ck 
下 
UU DQP 3 28? - SI +6) e S Solto + altok); 


(Gii) R2DPCA FU 2595 P, 的 均 方 误差 为 


$ N 
E-ri = SIL Y Yu 


DONA 3 x2. m m mnm m 
Neid E ee NSSS 9777 


m 


m n T 
ULTIO FON TLEN 9 3 Y Y edet otot), 


i1 j=1 121 k=1 
证 明 (i) Teh (ij) 7t 
tij = yes 
它 的 估计 为 


gs S S M 
i = i bx is" js? 


s=1 l=1 


(14) 


(15) 
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E(bi, b b by,) = 


B(S Yb, bt bibi) =E HED 


t 
ll 
pa 
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N-D Eyo! al, — ui) (o5, — uS) (o — ut) (od. — ui) 


n: E(o5, — ui)o5. — u5)( 2 (ofa =) ( x (0$ - 1) 


tips irs 
N -1y X x 
&O CIE pol, — ual, -u Y; (js i) È (Ofa = D) 
rem rm 
(N = 1)? L $ k k A u 
ume Elois = Ki)(Qjs 一 i 27 (os, 一 i5) ( » (os, — ui ut)) 
Een ts 
N —1y x 
(CD. pal, x i) (ot, 一 ue 5 (ai。 一 Hi i) ( D (o&, — Hi 3) 
sire irs 
N= D’ r l TA" e NI 
+a Elas - usas — a Y (ois 一 中) Y: (oi, — Hi ut) 
six zs 
N-1y x al 
OU Na V E(of, — uf) (os, — is >》 ois 一 n))( Y o. - ij) 
Ki" s 
1 N N N 
k k 
+E( > Qis 一 DS (9s - ub ))( Xe, - ub)( 3 oss - u5)) 
eire cre er ir 
nd 2 ny -— 
= DOS ONCE vae, QI UON D olot + olot + olto), Q6) 


: (ol, — ul) E m 内 ) 


3 一 1 


i2. ix: 
N 
N-1 1 
xi N (ad, EXE u$) P2 N »» (Qs o i5) 
s273 
N 
(5x - Qit 二 Hi ut) 
t1=1 
tst 
N 
N-1 1 
GC - ub - t E oa- ii 
tz 
_ (QN - UQN - 3) Vili 十 (N - 1} (N° - N? -2N +1) EN 
N2 ijij xxm 0i5;0ij 


QU -DCUV X SN S BNCESI Ies p oltott), 


N^4 Oii Ojj + Tij Tij 


(17) 
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因此 
" T m , T m , ^ —Na 2 m m 2 
El. 一 Te 人 = 3,9 EG - t) 2 >, EQ 十 N 359305) 
i=1 j=1 i-1 j=1 i=1 7 一 1 
= 二 Y e(3ohatal, «37 ta) 
i=1 j=1 l=1 k-1 s=1 m 1 


TEX (16) 和 (17) 代入 上 式 再 整理 即 证 。 
(ii) 的 证 明 与 (这 类 似 。 
3.2 样本 协 方差 阵 的 估计 准确 度 分 析 
下 面 比较 PCA 与 C2DPCA 样本 协 方差 阵 的 估计 准确 度 。 考 虑 光 与 六 的 均 方 误差 的 差 
A(E-f.) = EÈ- El - Elf. -Telt = = OU 3x x (V — Vds) 


1 1=1 k=1 


i—l]j 


E 3 2i m m n n 
二 


N SSN EN AEN =i z 
(Wo —— 2222222, (8e - ese). (18) 


看 一 个 例子 。 设 2 阶 随 机 和 矩阵 


满足 1 1 1 
£il = z372 = 372 一 4222 五 (zll) 一 0， E(x?) = 1, E(z14) => 0. 
经 计算 可 得 
AA AN 一 40N* + 32N? — 328N? + 456N — 192 
A C (en SN I E3287 NT N (ag 


车 取 ”n 05, WAN > 12 时 ， 式 (19) 小 于 零 。 这 个 例子 说 明 ，C2DPCA 样本 协 方差 阵 的 估 
计 准 确 度 未 必 一 定 优 于 PCA 样本 协 方差 阵 的 估计 准确 度 。 换 句 话 说 ，2DPCA 样本 协 方差 阵 的 
估计 准确 度 要 优 于 PCA 是 需要 条 件 的 。 

定理 4( 判 定 条 件 D 1) 车 


W EIE. — Tz < EIS — xl 
2) 若 


mn m 
2 XS 
i=1 i=} 
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JW EJt.-r.z«rzjE£-x(|2. F t > tz >- > ta ÆT, 的 特征 值 。 
证 明 将 式 (18) 中 的 1,k 互 换 ， 再 与 式 (18) 相 加 ， 得 


2 n n 


A x N ES 1 m m 
2A, o = a DDD VAS evt m 


i=1 j=1 l=1 k=1 


(N —1)( N3 + 2N? — 8N +6) Any lk „kl lk kl 
zi N6 25 2, 2 (oi i + 0i70 oij — 20 Tij oij) 
i=1 j=1 l=1 k=1 


Et 3 2. m m n n 
HOO EDD (tel toko- olot) 


i=l j-1i-—1k-1 


N*-3N* +N? -4N- 1 NN 
u—— — 99 99 logo ogo ooy) (20) 


iiza =d 地 一 由 ， 它 是 归 零 化 随机 图 像 矩 阵 4 — ECA) KI (4,7) zo, Ul 


lklk klkl llkk 2 2.2 
Viu + Viji 4 = E(x 1222.) F E(xvivj) — 2E(xüzjtikTjk) 


= E(züzjy zjzu)) 20, (21) 
0150 + 01501; — 201,01 = (E(zazjk)) + (E(zixzji))) — 2E(zazjk) E(zaz jt) 
= (E(zuzjk) 一 至 (zikzi)) > 0, (22) 
0405] + 05 0); — 20505 = E(zz) E(12,) + E(z,)E(z2)) 一 2E(zaza)E(vjij) 
2 2V E? (raza) E? (2j) — 2E (zuzik) E(Tji£jk) = 0, (23) 
即 式 (20) 的 前 3 项 非 负 ， 最 后 一 项 的 求 和 部 分 为 


n 


m m TL 
Hi Fk lk lk kl kl 
> > > > ogoi - odo - oot) 


i=1 j=1 1-1 k-1 


-(YY4-YYre)-Q4-Y1390. ea 


I= j—1 ty 


即 1) 成 立 ， 同 理 可 知 2), 也 成 立 。 
定理 5 (判定 条 件 TT) 若 
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证 明 HE3, D 5 D. 的 均 方 误差 的 差 为 


i=l j=1 l= 


i 3 2 m mn mn 
UN DON tad "T ol ie allg. kk) 


i=1 j=1 l=1 k=1 
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fr BernDBS 人 脸 库 和 JAFFE 人 脸 表 情 库 中 进行 人 脸 图 像 识别 实验 来 验证 上 述 分 析 。BernD 
BS 人 脸 库 由 30 人 组 成 ， 每 人 10 幅 ， 含 有 一 定 程 度 的 人 脸 姿态 、 表 情 和 饰物 的 变化 ， 每 幅 图 
1879 128 x 128 的 矩阵 。 每 人 随机 取 2 幅 共 60 幅 用 于 训练 ， 余 下 用 于 测试 。JAFFE 人 脸 表情 库 
由 10 名 目 本 女性 组 成 ， 含 有 生气 、 厌 恶 、 害 怕 上、 高 兴 、 中 性 、 悲 伤 和 惊奇 7 种 表情 ， 每 人 每 种 
表情 若干 幅 ， 共 213 幅 表情 图 像 。 手 工 标定 两 眼中 心 位 置 ， 并 进行 校正 和 裁剪 ， 将 表情 图 像 处 
理 成 128 x 128 的 图 像 矩阵 ， 每 人 每 种 表情 取 工 幅 共 70 幅 作为 训练 样本 ， 余 下 的 143 幅 作为 测 
试 样本 。 每 个 库 上 的 实验 都 进行 10 轮 ， 取 10 轮 结果 的 平均 作为 最 后 的 结果 。 

分 别 采 用 PCA 和 2DPCA 进行 人 脸 图 像 识 别 ， 采 用 最 近邻 分 类 器 进行 分 类 ，PCA 采用 2- 范 
数 作为 距离 度量 ，2DPCA 采用 Frobenius 范 数 作为 距离 度量 。 在 实施 PCA 与 2DPCA tt, Pr 
取 投 影 向 量 的 个 数 k 由 投影 特征 值 比例 来 确定 ， 即 可 用 下 式 来 确定 


k 
k=argmin{ 5 di/d>e), (25) 


其 中 由 可 以 是 Xi， 先 或 如 ,qd 为 从 的 特征 值 之 和 ， 其 实 也 是 信和 全; 的 特征 值 之 和 ，c € [0, 1]. 
表 1 给 出 了 PCA 样本 协 方差 阵 与 2DPCA 样本 协 方差 阵 的 特征 值 平方 和 和 ， 图 1 与 图 2 分 别 给 
出 了 BernDBS 库 与 JAFFE 库 的 实验 结果 。 对 照 表 1 与 图 1，2 可 见 ， 实 验 结果 与 定理 4、 定 
理 5 的 结论 相符 。 
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图 1: BernDBS & P PCA 5 2DPCA 的 比较 图 2: JAFFE &'* PCA 5 2DPCA 的 比较 
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表 1: PCA 与 2DPCA 的 样本 协 方差 阵 特 征 值 平方 和 的 比较 


特征 值 平方 和 PCA C2DPCA R2DPCA 


16384 , 128 128 
BernDBS # ÀZ = 5.9666 x 1014 $552 = 8.7427 x 10" y; 8 = 1.2945 x 1015 
EFI i=l =l 


i= 
16384 


n 128 
JAFFE 库 和 2 — 6.9073 x 10? 412 = 1.0189 x 101^. | © 2 = 1.1340 x 104 
i ini 


i=1 


5 ”结束 语 


定理 4 给 出 了 2DPCA 在 图 像 特征 提取 中 优 于 PCA 的 判定 条 件 ， 是 一 个 充分 条 件 ; 定 
理 5 给 出 了 C2DPCA 与 R2DPCA 在 图 像 特 征 提取 中 熟 优 熟 劣 的 判定 条 件 ， 是 一 个 充 要 条 件 ; 
所 有 这 些 条 件 都 只 依赖 于 各 自 协 方差 阵 的 特征 值 平方 和 的 大 小 。 
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Criteria for 2DPCA Superior to PCA in Image Feature Extraction 
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Abstract: PCA and 2DPCAÀ are two typical image feature extraction methods. The wellness of image 
features extracted by PCA and 2DPCA can be judged by their reconstruction errors. Our analysis 
shows that the reconstruction errors of PCA and 2DPCA are equal in theory, but depend on the 
estimation accuracy of their covariance matrices in practice. This paper gives the accuracy expressions 
for estimations of covariance matrices for PCA and 2DPCA, respectively, based on the criterion of 
mean square error. Further, it is shown that the criterion for 2DPCA to be superior to PCA is that 
the square sum of eigenvalues of the covariance matrix of 2DPCA is larger than that of PCA. This 
paper also points out that the same criterion can be used to judge the superiority between the column 
2DPCA and row 2DPCA. Experiments on the face database and facial expression database verify the 
correctness of the decision conditions. 

Keywords: PCA; 2DPCA; image feature; covariance matrix; reconstruction error; mean square error 


